Universality and backreaction in a general-relativistic accretion of steady fluids 
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The spherically symmetric steady accretion of polytropic perfect fluids onto a black hole is the simplest flow 
model that can demonstrate the effects of backreaction. The analytic and numerical investigation reveals that 
backreaction keeps intact most of the characteristics of the sonic point. For any such system, with the free 
parameter being the relative abundance of the fluid, the mass accretion rate achieves maximal value when the 
mass of the fluid is 1/3 of the total mass. Fixing the total mass of the system, one observes the existence of two 
weakly accreting regimes, one over-abundant and the other poor in fluid content. 
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The spherical steady accretion of gas onto a gravitational 
center has been investigated in the Newtonian context by 
Bondi 1 1 ] and in the Schwarzschild space-time by Michel 1 2 ] , 
Shapiro and Teukolsky |3], and others |4]. The Newtonian 
critical accretion flow was shown to be stable by Balazs and 
others |5], while the stability of the steady flow in the fixed 
Schwarzschild background has been analyzed by Moncrief 
1 6]. The general-relativistic description including backreac- 
tion has been formulated in [8J]. There are two reasons for 
inspecting steady flows in the full generality, when infalling 
gas modifies the space-time geometry, that is including back- 
reaction 1 7]. First, one can see the effects of backreaction in 
a simple but nontrivial accretion model. Second, gravitational 
collapse of the fluid which starts as a flow dominated by a 
steady accretion, can be better controlled, and that would al- 
low one to have an insight into formation of gravitational sin- 
gularities. This paper focuses on the issue of backreaction. 
We show that the backreaction strongly influences the mass 
accretion rate. 

We will consider a spherically symmetric compact ball of 
a fluid falling onto a non-rotating black hole. The black hole 
provides the simplest choice for the central object since one 
can neglect the occurrence of shock waves. The fluid is re- 
garded to be steady in a sense defined later and by a black 
hole we understand the existence of an apparent horizon. For 
the detailed derivation of equations the reader should consult 
0. Here we give only a brief description, focusing attention 
on the physical assumptions. We will use comoving coordi- 
nates 



ds 2 = -N 2 dt 2 + adr 2 + R 2 (d0 2 + sin 2 6d(p 2 ) , 



(1) 



where the lapse N, a and the areal radius R are functions 
of a coordinate radius r and an asymptotic time variable t. 
The nonzero components of the extrinsic curvature Kn of the 
t = const slices read K r r = ^d,a, K* = K e g = |§ = 

^(trA" - K r r ). Here trK = j^d, In ( ^R 2 }. The mean curvature 
of two-spheres of constant radius r, embedded in a Cauchy 
hypersurface is p - j^j= . The energy-momentum tensor of 
the perfect fluid reads T^ v = (p + g)^^ + pg^ v , (u^uf 1 = -1) 
where u M denotes the four-velocity of the fluid, p is the pres- 
sure and q the energy density in the comoving frame. 



The areal velocity U of a comoving particle designated by 
coordinates (r, f) is given by U(r, t) = ^ = j(trK-K r r ). From 
the Einstein constraint equations one has 



(2) 



and d R (R 2 U) - R 2 trK = 0. We use here and in what follows 
the relation d r = y/a(pR/2)d R in order to eliminate the differ- 
entiation with respect the comoving radius r. The quasilocal 
mass m(R) obeys the equation dRin(R) = 4jtR 2 q. The mass 
accretion rate is m = (d t - (d t R)d R )m(R) - -4nNR 2 U(g + p). 
A more familiar form of that is m = -4nR 2 nU, where n is 
the baryonic density (see below). A standard condition for the 
steady collapsing fluid is that all its characteristics are con- 
stant at a fixed R (see |9]). In analytical terms (5,X|« =const — 
(d t - (d,R)d R )X = 0, where X = o. U.a. . . ffoll. The Einstein 
evolution equation d,U = p 2 R 2 d R N/4 - m(R)N/R 2 - 4nNRp 
and the energy conservation equation d,g = -NtrK{g + p) be- 
come ordinary differential equations with respect R, assuming 
steady flow. Selfgravitating steady fluids still allow for the 
change in time of some geometric quantities, like the mean 
curvature p or the area of the black hole. It is easy to show 
that for a steady flow d R m — 1 8]. 

By a black hole is meant a region within an apparent hori- 
zon to the future, i.e., a region enclosed by an outermost 
sphere S a on which the optical scalar 0+ = ^- + U vanishes 
ill . (The other condition, that 6L = Si _ jj > o for all 
spheres outside Sa, is satisfied trivially for a steady accre- 
tion.) That means that on Sa the ratio 2m BH /R A H becomes 
1 , where Rah is the areal radius of the apparent horizon and 
niBH = rn(R.AH) is the mass of the black hole. We shall spe- 
cialize to polytropic perfect fluids p = Kg r , with F being a 
constant, (1 < T < 5/3). 

Furthermore, assume the radius R^ of the ball of fluid and 



boundary data are such that \U X \ <sc 



<s a„. These 



boundary conditions are needed in order to glue the steady 
fluid with the external Schwarzschild geometry (see a discus- 
sion later). The momentum conservation equation V^T^ = 
— in comoving coordinates 



Nd R p + (p + g)d R N = 



(3) 



2 



can be integrated, yielding (with NiR^) =1) 

r + al. 



-r + 



(4) 



Here k = ^X. @ can be regarded as the general-relativistic 
version of the Bernoulli equation. Another useful relation that 
follows from Q is N = jj^ . Here appears the baryonic num- 
ber density n — C exp ds s +p is ) ■ The baryonic number den- 
sity works, for perfect barotropes, as an integration factor; if 
T^y is conserved, then V„ (nu^) = 0. This representation of 
the lapse function shows that the lapse N is constant along an 
orbit of the constant areal radius R. 

In the general-relativistic case the sonic point is defined 
as a location where the length of the spatial velocity vec- 
tor \U\ - \U\/(pR/2) equals a. Therefore at a sonic point 
\U\ - \pRa. In the Newtonian limit this coincides with the 
standard requirement \U\ = a. In the following we will de- 
note by the asterisk all values referring to the sonic points, 
e.g. a*, U», etc. The four characteristics, a„ {/», m*/R, and 
c« are related H, al (l - ^ + c*) = Ul = ^ + c„ where 
c* = 2nRlp*. The infall velocity U reads 



U = U *f> 



1 + -T 

a- 



i/(r-i) 



(5) 



Here t/» is the negative square root. From the relation be- 
tween the pressure and the energy density, one obtains, using 
equation @ 



2/(T-l) 



= g a 



r 



4 + P 



N K 



(6) 



where the constant Q m is equal to the mass density of a col- 
lapsing fluid at the boundary The steady fluid is described 
by equations (|4]-|6j. They constitute an integro-algebraic sys- 
tem of equations, with a bifurcation point at the sonic point, 
where two branches (identified as accretion or wind) do cross; 
that is a well known feature of that problem, present also in 
models with a test fluid 111 - 15], 1 8]. That requires some cau- 
tion in doing numerics and a careful selection of the solution 
branch. Notice that these equations are expressed exclusively 
in terms of quantities that are steady in the sense of the former 
definition. 

This collapsing compact cloud can be connected to the 
Schwarzschild exterior. The ideal model of a steady flow 
would be that engineered by a "daemon" who just fixes bound- 
ary data conditions; a similar picture might well be valid in 
some binaries. Alternatively, one can imagine an initial con- 
figuration with a steady fluid annulus around the black hole 
and an external transition zone that is not steady (since both 
pressure and the mass density must be made to vanish there) 
and that can expand outward and inward with the speed of 
sound, perturbing the external layers of the steady interior. 
The boundary can be connected at via an almost mass- 
less transition zone, to the external Schwarzschild geometry; 



the standard junction conditions — continuity of the metric 
and of the transversal extrinsic curvature, and the condition 
Noo — 1 — require now the appropriate choice of the slicing 
of the exterior region. We verified in a number of numerical 
calculations that for short evolution times (f <K /?oo/flco) and 
under the former boundary conditions the bulk of the fluid is 
steady. (Notice, that this observation indicates that such ini- 
tial data can be useful in order to model the gravitational col- 
lapse.) The asymptotic mass m of the spacetime would be well 
approximated by moo, the sum of the masses of the black hole 
and the quasi-stationary gas, for a suitably chosen transition 
zone. 

In numerical calculations it will be convenient to represent 
the mass m(R) in the form 



m(R) 



= m — 4jt I 

Jr 



drr 2 g. 



(7) 



Notice that m(oo) = * m. Assuming steady flow, one can 
obtain following expression for the lapse N |H 



pR 

N = .1 — P(R), 



p(Rco)R. 
P(R) = exp 



,_ , ds 

(p + e)-r 

p L s 



(8) 



The full description of the steady accretion that will be dis- 
cussed later is given by equations ©, 0, 0, l|3 and ©; this 
is equivalent to the former system. 

One of two main results of this paper is realization that sig- 
nificant information about the full system with backreaction 
can be obtained through the investigation of steady flows with 
the backreaction being ignored. 

Namely,let be givenF e (1,5/3], asymptotic data and 
cioc and let \Uco\ <s <k a^. Assume that the sonic point 
is located outside the apparent horizon, i.e., XJ\ < 0.25. 
Then the sonic point parameters a», U\ and m t /R t in the 
above model with backreaction are the same as in the test 
fluid accretion with the same asymptotic data 1 11]. 

Sketch of the reasoning. We will apply the bootstrap-type 
argument. Let us assume temporarily that C/J < 0.1, that is 
R t > 5m«. This implies pR/2 > VOJ, a 2 < al < 1/7 and 

\6tt ^{p + e )jrds < 80^2^ for R > R t . 

If m - m(R*) < (49/80)m(^»), then the integrand present in 
the exponent of the function P(r) is smaller than m*/R*. That 

in turn means that > 

^1 _ 4*i + ul and 



at < -r + 



r + al 



(9) 



Solutions of this inequality are bounded from above by the 
solution of the equation 



a 2 = -r + 



r + al 



(10) 



3 



where to = 2m*, at the point where a 2 / (l + 35 2 ) = m/(2R*) 
I . Note that R* > 2m; this allows us to use the estimate of 
Theorem 2 of 1 8]. Thus one has 



ing fluid can be expressed as below (see Eq. (6.1) in |8J) 



at < a < 



2al 



5 -3T + 



3a 2 (r-l) 2 (9r-7) ' 
4r(l+3a 2 ) 



(ID 



Inserting this estimate into 0, one can bound the matter den- 
sity, g(R) < goo (l + C(T, a,)/ (5 - 3r + a 2 )). Here the con- 
stant C is of the order of 10. This in turn allows one to replace 
the bound on the exponent of /3 by a much better one of the or- 
der (m*/R*) (Rt/Rca); this is clearly negligible in comparison 
to m*/R* (save for a small region around r = 5/3) and one can 
repeat the previous argument, but this time assuming /3 - 1 . 
This would lead to the improved estimate of a and g outside 
R = 5R*. Further refinement — the extension into bigger val- 
ues of U 2 and of m*/m, and to all F's in the interval (1,5/3] 
— would be done by splitting the integral of ft into two inte- 
grals, °° . . . = J R *... + § 5R °* . . . and separately estimating 
each term. 

Using the preceding estimates, one gets 



C = -2^ < Cal-^ « -. 



(12) 



That implies that both c* and the exponent of ft can be put to 
zero in all relations between characteristics of the sonic point. 
Inserting this information into the Bernoulli equation @ one 
obtains that at a sonic point 1 8] 



i+y(3r-i) = 3(<£ + r) 



y 



(13) 



where y = 1 — 3m.J(2R t ). Coefficients of this algebraic equa- 
tion do not depend on the asymptotic mass density and there- 
fore y is ^-independent. The sonic mass m* and the sonic 
radius R* clearly depend on g^ but their ratio is constant. In 
fact, m*/R* must be exactly the same as in the case when the 
backreaction can be neglected, that is when the mass of the 
fluid outside the black hole is small in comparison to the total 
mass. The same conclusion holds true also for other parame- 
ters of the sonic point, the fluid velocity U* and the speed of 
sound a*. In conclusion: a 2 , U 2 and m*/R* can be inferred 
from a suitable steady flow with a test fluid. 

It follows from the above discussion that the mass den- 
sity changes moderately outside the sonic point. In fact, for 
R € [R*, 10R*] g(R) < C(R)g cx> with C(R*) being of the or- 
der of 10 and C(10/?«) being slightly larger than 1. Assum- 
ing that R m » 10R t one can show (taking into account the 
fact, that the region (10/?*, /?«>) has a much larger volume than 
(R*, 10/?*)) that under the preceding assumptions the mass of 
outer layers of the fluid m - m* = 4tt J r drr g = yg^ with a 
constant y determined in practice only by R m . Alternatively, 
m*/m = 1 — gca(y/m); the ratio m*/m is a linearly decreasing 
function of . 

The rate of the mass accretion m within the steadily accret- 



Rl 



m - nm^oo—z —r- 



(5-3D 

-,2 \ 2(r-n 



a 2 \ 1 + 3al 

1 + 



(14) 



The whole dependence on gco is contained in the factor m^g x . 
This achieves (fixing m) a maximum at m* = 2m/3 and tends 
to zero at both ends: i) m* — > m (when the density g tends 
to zero) and ii) m,/m — > (when the mass of black hole is 
negligible in comparison to the mass of the fluid). In this place 
we invoked the assumptions concerning boundary conditions. 

Thus, one arrives at our second main result: Amongst 
steadily accreting systems of the same and m those 

will be most efficient for which m, = 2m/3. The factor 2/3 
is universal — independent of the parameters Ra,, T and a^. 

Since the mass of the fluid m/ is close to m - m» and m = 
ntf + niBH, the above means that the maximum of the mass 
accretion takes place when the mass of the fluid is a half of 
the mass of the black hole. We would like to point that in 
H41 the mass to, is only approximately constant, and hence to 
is not strictly constant in time. The fulfilment of tot «: m„, 
where r is the evolution time, can be understood as yet another 
criterion for the validity of the steady flow approximation. 

In the remaining part of the paper the asymptotic mass to 
has been normalized to 1 . In order to test the above results, 
we performed a series of numerical calculations, according to 
the following scenario, i) Fix a value of the adiabatic index F 
and the asymptotic speed of sound a\\ change in a systematic 
way the asymptotic mass density, ii) For the same F choose 
new a 2 ^ and change the asymptotic mass density, iii) Choose 
new T and repeat the preceding procedure. Below we shall 
describe only a small subset of the numerical evidence. 



TABLE I: Characteristics of the sonic point a 2 t , |t/»|, R, for the T = 
4/3 poly trope and a 2 m = 0. 1. The last column shows the areal radius 
of the apparent horizon. The mass accretion rate (second column) 
achieves maximum at m BH ss 0.663. 

nif/m m R, (J, a~ Rah 



4.186 x 10 
4.186 x 10 
0.041 
0.083 
0.167 
0.251 
0.333 
0.418 
0.628 
0.837 
0.879 
0.962 
0.983 



4.156 x 
4.156 x 
3.815 x 
6.978 x 
1.152 x 
1.397 x 
1.469 x 
1.403 x 
8.612 x 
2.187 x 
1.264 x 
1.277 x 
2.382 x 



io- 48 
io- 33 
io-' 8 
io-' 8 
io-' 7 
io-' 7 
io- 17 
io- 17 
io-' 8 
io- 18 
io- 18 

10 19 

1Q -20 



4.290 
4.290 
4.110 
3.931 
3.571 
3.212 
2.859 
2.493 
1.594 
0.696 
0.516 
0.156 
0.067 



0.34137 
0.34137 
0.341369 
0.341371 
0.341371 
0.341371 
0.341371 
0.341399 
0.341369 
0.341394 
0.341399 
0.34137 
0.341399 



0.179184 
0.179184 
0.179182 
0.179183 
0.179182 
0.179183 
0.179183 
0.179184 
0.179182 
0.179185 
0.179184 
0.179187 
0.179185 



1.993 
1.993 
1.906 
1.821 
1.656 
1.492 
1.326 
1.156 
0.740 
0.322 
0.239 
0.072 
0.031 



It appears that the characteristics a,, U 2 , m t /R t of the sonic 
point practically do not depend, for a given F and a 2 ^, on 
the energy density (exemplary results are shown in Table 



4 



1), suggesting a precision of 10~ 5 . The more reliable mea- 
sure of the numerical error is the difference between the mass 
m = 1 (assumed in the equations) and the mass found numer- 
ically. Our results (see Table 1) suggest that numerical error 
is smaller than 0.5 %. The quantity c, is negligible in com- 
parison to other sonic point parameters. The mass m, behaves 
like m, = 1 - yg^ (Fig. 1). And finally, the extremum of m 
is achieved when m, ~ 2/3 (Fig. 2). Notably, numerical re- 
sults extend the regime in which the universality is observed 
into fluids characterized by the speed of sound a x exceeding 
1. This is remarkable, since that implies R* < Rah and the 
theoretical analysis presented above would not be valid. 
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FIG. 1: The dependence of ~ 1 - m, on the asymptotic mass 
density. 
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FIG. 2: The dependence of m on the ratio (m - m,) /m. 

Formula (I14> implies that the mass accretion rate in the case 
with backreaction (when m„ < m) can be only smaller from 
the mass accretion rate without backreaction (when m, = m), 
with the same asymptotic parameters p^,, a M , m. Numerical 
results confirm the preceding argument that y depends essen- 
tially only on the size of support. 

In the case of polytropes with p - Kn Y the Bernoulli equa- 



tion reads N\T — 1 - a^j = T- 1 -a 2 . The analytic arguments 
of the type presented above can be applied and conclusions are 
the same. Similar results should be true for Newtonian mas- 
sive fluids for F significantly smaller than 5/3. 

In conclusion, in the simple model of accretion with back- 
reaction considered here, one can get all parameters describ- 
ing the sonic point, with the exception of its location R* and 
mass m* simply from a related polytropic model with the test 
fluid having the same index F and the same asymptotic speed 
of sound Oco. The main result is that the mass accretion rate 
m achieves maximum at mf/mgH ~ 1 /2; therefore there exist 
two different regimes, mf/niBH ^ 1 and mf/mgH » F with 
low accretion. In contrast to that, in the test fluid approxima- 
tion the quantity m grows with (that is with mf/niBH)- The 
mass accretion rate m sets the upper limit for the luminosity 
of the system. Hence this simple analysis of accretion sug- 
gests the existence of two weakly luminous regimes: one rich 
in fluid with mf/msn s> 1 and the other with a small amount 
of fluid, nif/mBH ^ F We noticed also that the normalized 
efficiency of accretion (defined as m/niBH) goes to zero. 

Results of this paper show the importance of the backre- 
action. We believe that the qualitative features demonstrated 
by the spherically symmetric model - the dependence of the 
mass accretion on mf/niBH and the existence of two weakly 
accreting regimes - will appear in the descriptions of accret- 
ing systems with angular momentum. 
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